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Ql,; ABSTRACT 

The production of runaway massive binaries offers key insights into the evolution of close 
+3 , binary stars and open clusters. The stars HD 14633 and HD 15137 are rare examples of such 

runaway systems, and in this work we investigate the mechanism by which they were ejected 
from their parent open cluster, NGC 654. We discuss observational characteristics that can 
be used to distinguish supernova ejected systems from those ejected by dynamical interactions, 
and we present the results of a new radio pulsar search of these systems as well as estimates of 
in ' their predicted X-ray flux assuming that each binary contains a compact object. Since neither 

pulsars nor X-ray emission are observed in these systems, we cannot conclude that these binaries 
contain compact companions. We also consider whether they may have been ejected by dynamical 
interactions in the dense environment where they formed, and our simulations of four-body 
interactions suggest that a dynamical origin is possible but unlikely. We recommend further X- 
ray observations that will conclusively identify whether HD 14633 or HD 15137 contain neutron 
stars. 

Subject headings: binaries: spectroscopic, stars: early-type, stars: kinematics, stars: winds, stars: indi- 
vidual (HD 14633, HD 15137), open clusters and associations: individual (NGC 654), pulsars: general, 
X-rays: binaries 

1. Introduction some are observed at high galactic latitudes and 

with large peculiar space velocities. The estimated 
While most O- and B-type stars are believed fract ion of O stars that are runaways range s from 

to form in open clusters and stellar associations, 7% (jConti. Leen. fc Lorrej fl977h to 49% jStone 



1979). These stars were likely ejected from the 



1NSF Astronomy and Astrophysics Postdoctoral Fellow clusters of their birth, and there are two accepted 
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mechanisms to explain the origin of their run- 
away velocities. Close multi-body interactions 
in a dense cluster environment may cause one 
or more stars to be scattere d out of the region 
(jPoveda, Ruiz, fc Allenlll967l) . On the other hand, 
a supernova explosion within a close binary may 
eject the secondary star due to the conservation of 
momentum of mass l ost in the exp losion, includ- 
ing supernova kicks (|Zwickvlll957t 
Sutantvo|[l97i) . 
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The production of runaway O- and B-type bina- 
ries is expected to be rare, but these systems can 
offer key insights into the evolution of close binary 
stars and open clusters. Distinguishing between 
the dynamical ejection and binary supernova sce- 
narios for an isolated runaway star can be nearly 
impossible unless the star contains CNO-processed 
gas received from a past episode of mass trans- 
fer ijHoogerwerf. de Bruijne. fc de Zeeuwl 120011 ). 
However, the ejection mechanism responsible for 
producing runaway binaries can be distinguished 
more easily from a variety of observable proper- 
ties. 

Leonard fc Duncan ( 1990l ) use TV-body simula- 



tions of open clusters to show that a binary fre- 
quency of about 10% is expected from dynamical 
ejections. Most of their simulated binaries have 
eccentric orbits, 0.4 < e < 0.7 , and high mass ra- 
tios, q = MijM\ > 0.5. While lLeonard fc Duncan 
(|1990h do not discuss the resulting orbital peri- 
ods in their simulations, we do not expect long 
period (P > 1000 days) binaries to form via dy- 
namical interactions since their low binding ener- 
gies would likel y cause them to be io nized. On 
the other hand, iPorteg ies Zwartl (|2000h predicts a 
higher binary fraction of 20—40% among runaways 
which are ejected by supernovae in close binaries. 
The resulting neutron star may remain bound to 
the secondary if not enough mass is lost during the 
explosion or if the explosion is not symmetric (pro- 
ducing a kick velocity). In such cases, the present- 
day orbital period and eccentricity depend on the 
mass e jected during the supernova and the kick ve- 
locity dNelemans. Tauris fc van den Heuvell 19991 : 
Brandt fc Podsiadlowskilll99a ). but the resulting 
period distribution is e xpected to be 3—1000 days 
(|Portegies Zwartll2000l ). 

Distinguishing between the dynamical or super- 
nova ejection mechanisms among runaway spec- 
troscopic binaries can be a difficult task since 



both are expected to produce eccentric, relatively 
short period binaries. Only dynamical interac- 
tions are expected to produce double-lined spec- 
troscopic binaries (SB2s), while a single- lined sys- 
tem (SB1) may be formed either way. Obtain- 
ing very high signal-to-noise (S/N) spectra or per- 
forming doppler tomography can, in some cases, 
identify optical companions in such systems. How- 
ever, it can be nearly impossible to detect a cool, 
low mass, optical companion with an O- or B-type 
primary. 

Neutron star companions can be difficult to 
identify as well. A few runaway O stars are pre- 
dicte d to be associated with pulsars (IPorteeies Zwartl 
2000h . but so far none have been f ound (jPhilp et ahl 
19961 ISaver. Nice, fc Kaspi 1996f ) . X-ray emission 
may also be used to pinpoint neutron star com- 
panions since a neutron star may accrete mass 
from its companion via a Roche lobe overflow 
stream or by Bondi-Hoyle accretio n from the stel- 
lar w i nd of the luminous primary dBondi fc Hovld 
1944 iDavidson fc Ostrikerl 1 19731 lKaperl 1T998h 



The systems experiencing Roche lobe overflow 
usually have large X-ray luminosities and strik- 
ing optical emission lines. On the other hand, 
weaker X-ray emission will result if stellar winds 
are the accretion source. According to the wind 
accretion model, the observed X-ray luminos- 
ity depends on the system separation, the rel- 
ative wind and companion velocities, the stel- 
lar mass-loss rate, and the mass of the acc retor 
( Lamers. van den Heuvel, fc Petterson 19761 ). If 
the stellar wind is too weak, or the system sepa- 
ration too large, the X-ray luminosity would di- 
minish below the detection limit of today's X-ray 
surveys, hence a "quiet" massive X-ray binary 
(MXRB). 

We have recently completed a radial velocity 
study of 12 field and runaway O- and B-type stars 
to identif y spectroscopic binari es with high space 
velocities ([McSwain et al.ll2007l ) . In that work, we 
improved the orbital elements of five SBls, but 
among these only two are definitely runaways: HD 
14633 and HD 15137. HD 14633 is classified as an 
ON8.5 V star, and HD 15137 is somewhat more 



evolve d as an 09.5 III(n) star. iBovajian et al 



(2005) 1'ound that these runaway SB1 systems were 
likely ejected in separate events from the open 
cluster NGC 654, but the travel times since their 
ejection (14 Myr for HD 14633 and 10 Myr for 
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HD 15137) present an unusual paradox since they 
are longer than the expected lifetimes of O-type 
stars. They attempted doppler tomographic sep- 
aration for both systems, but this search for an 
optical companion was unsuccessful. HD 14633 is 
a nitrogen-strong star (|Walbornlll972l ). and Boya- 
jian et al. argue that its nitrogen enrichment and 
fast rotation resulted from a mass transfer episode 
prior to a supernova in the close binary system. 
They propose that both stars are quiet MXRB 
candidates that harbor neutron star companions 
but are not known X-ray emitters. 

To investigate this claim, we present here the 
results of a pulsar search in these two systems. 
We also compare the observed upper limits of their 
X-ray fluxes to model predictions of X-rays gener- 
ated by the accretion of stellar winds onto a com- 
pact companion. Since neither pulsars nor X-ray 
emission are detected, we are not certain that ei- 
ther system contains the product of a supernova. 
Therefore we also investigate the possibility that 
dynamical interactions may have been responsible 
for their ejection from NGC 654. 

2. Constraints on the Secondary Masses 

We recently presented updated orbital elements 
and mass functions, f(m), for the SB1 system s 

~' 2007f ). 



HD 14633 and HD 15137 (McSwain et al 



and these parameters are summarized in Table Q] 
While f(m) does not strictly define the component 
masses, it does allow constraints on the primary 
star's mass, Mi, the companion mass, M2, and 
the system inclination, i. We plot the relation- 
ships between these properties in Figures [T] and 
[2] using a grid of values for M\. In addition, the 
vertical dashed line in both figures indicates the 
expect ed mass of each star based on the calibra- 
tion of iMartins. Schaerer. fc Hillierl (|2005T ). Both 
binaries have very low f(m) that suggest low mass 
companions (1— 3M©). Without directly detecting 
the companions, we must constrain their masses 
using another method. 

We used the statistical method o f lMazeh fc Goldberd 
(|1992l ) to determine the most probable q, and this 
relationship is plotted as a solid line in the top 
panels of Figures Q] and O Note that their tech- 
nique results in a somewhat lower probable value 
for i than the expectation value (i) ~ 60° since 
it accounts for observational selection effects. HD 



14633 and HD 15137 are particularly interesting 
since they probably have q < 0.3, and such low 
q are unknown among O star binaries (exclud- 
ing X-ray binaries). Their allowed masses do not 
exclude neutron star companions, and below we 
investigate the possibility of radio pulsars or a de- 
tectable X-ray signature of a compact companion. 

3. Pulsar Search Results 

HD 14633 was targeted in a pre vious search 
for pu lsars in runaway OB stars by et al. 

( 19961 ). Although they did not find any pulsars, we 
chose to perform a new, more sensitive search to 
inve stigate the companions of both runaway bina 



Portegies Zwart (2000) shows that detection 



ries. 

of pulsars in O star binaries is hindered by the ab- 
sorption of the radio emission by the stellar winds. 
This absorption will be greater at periastron than 
apastron since the winds are densest close to the 
O st ar, and our new ephe merides for these bina- 
ries (jMcSwain et al.ll2007l ) allowed us to schedule 
the observations as close as possible to the max- 
imum separation. (If these systems have high i, 
true if a neutron star is present, the times of infe- 
rior conjunctions may provide the lowest column 
depth instead. These times correspond to the or- 
bital phases 0.83 for HD 14633 and 0.70 for HD 
15137. Although our observations were not timed 
to occur at inferior conjunction, they are well sepa- 
rated from superior conjunction when the neutron 
star could possibly be eclipsed.) We also chose to 
use higher frequencies than the previous investiga- 
tions because they are more likely to detect radio 
pulses dispersed by the stellar winds. 

We observed HD 14633 and HD 15137 using 
the National Radio Astronomy Observatory's 100- 
m Green Bank T elescope (GBT) and the Pulsar 



Spigot back-end (jKaplan et alj 120051 ). We used 



the Spigot in mode 2 with the S-band receiver to 
obtain 600 MHz (1650-2250 MHz) of bandwidth, 
which the Spigot summed and synthesized into 768 
0.78125-MHz frequency channels every 81.92 /j,s. 
The Spigot in mode 42, with the 820 MHz re- 
ceiver, provided 50 MHz of bandwidth (795-845 
MHz), summed into 1024 0.0488-MHz frequency 
channels with the same time resolution. Each tar- 
get was observed once with each receiver for 6300 
s. For each observation, the time at midexposure 
and the corresponding orbital phase are listed in 
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Table [2j Although we requested observations at 
the times of apastron (orbital phase 0.5), telescope 
scheduling constraints limited us to slightly earlier 
times during each orbit. 

The data we re reduced usin g the PRESTO soft- 
ware package (| Ransom! 120011) . We first searched 
the raw data for radio frequency interference in 
both the time and frequency domains and applied 
interference masks to remove these sources. Each 
of our targets is expected to have a dispersion mea- 
sure (D M) ~ 50 - 75 pc cm" 3 using the NE2001 
model ( Cordes fc Laziol 2002 ) w ith the distances 
in Table Q] ( McSwain et al.1 2007 ). so we searched 
the observations by dedispersing the raw data into 
separate time series with DMs ranging from to 
120 pc cm~ 3 and spaced by 0.5 pc cm -3 . We then 
Fourier-transformed each time series and searched 
them by using Fourier-domain acceleration search 
techniques in order to maintain sensitivity to both 
pulsars in binary orbits as well as isolated pul- 
sars. For the acceleration search, this technique 
involves a linear approximation of the binary's 
orbital acceleration (valid since the observations 
span a small fraction of the orbit) to look for a 
derivative in pulsar frequency. The flux density 
sensitivity, Sn m , and luminosity, Li; m , limits of 
each observation are listed in Table [2] The typi- 
cal pulsar has a luminosity L ^ 0.1 mJy kpc 2 at 



use 



1400 MHz (jLorimer et al.N2006h . and our L lim do 



go fainter than this value for the S-band observa- 
tions. 

While we did not detect any pulsars in these 
systems, we cannot rule out their existence since 
a pulsar's beam may not cr oss our line of sight. 
Tauris fc Manchester! (|l998h find that the ob- 
served beaming fraction depends on the rotational 
period, P, since the beam radius varies as P _1//2 . 
P, in turn, depends on the pulsar's age, t, since 
magnetic dipole radiation causes the pulsar to spin 
down over time. Without detecting a pulsar we 
cannot measure P directly, but if the runaways 
were ejected by supernovae, their dynamical ages 
correspond to the ages of the resulting neutron 
star. HD 14633 was ejected from the open cluster 
NGC 654 about 14 Myr ago, and H D 15137 left the 



same cluster about 10 Myr ago (jBovajian et al 
2005). If we assume that any pulsar was born 



spinning much faster than its current rate, we can 



P 



8tt 2 P 6 
3c 3 / 



B 2 t, 



(1) 



where R is the neutron star radius, c is the speed 
of light, / is the moment of inertia, B is the 
strength of the mag netic field, and t is the pu l- 
sar age in seconds ( Tauris fc Manchester 1998h . 
Using typical neutron star values of R = 10 km , 
/ = 10 45 g cm 2 , and B = 10 12 G (lLipunovlll992h . 
we estimate P » 1 s for HD 14633 and HD 15137. 
Lori mer et al. (|2006l) find a beaming fraction of 
~ 20% for such pulsars, making it unlikely that 
we will detect the beam. 

Clearly, our GBT observations are not sufficient 
to characterize the companions in cither system. 
Below we investigate their other properties to de- 
termine if they should display an observable X-ray 
luminosity. 

4. Predicted X-ray Luminosities 
4.1. Wind Accretion Model 

HD 14633 and HD 15137 are SBls with high ec- 
centricity, but neither are known X-ray sources. If 
they do contain compact companions, should these 
stars have a detectable X-ray emission? We pre- 
dict the X-ray flux of each quiet MX RB candidate 



using the wind accretion model of lLamers et al 
(197a), modified for eccentric orbits. This method 
uses the Bondi-Hoyle accretion rate 



S a 



Ana 2 



(2) 



(La mers et al. I ll976h . which depends on the stellar 
mass loss rate, M, the system's major axis, a, and 
the accretion radius of the companion, 



Ra = 



2GM 2 



(3) 



Here, M 2 is the mass of the accretor and V re i is the 
relative velocity of the stellar wind and the com- 
panion's orbit, V W i n d and V or b respectively. V or b 
is determined from a and the orbital eccentricity, 
e. The predicted X-ray luminosity is 

L x oc {S a M 2 (4) 

( Lamers et al. 19761) . The efficiency, £, of convert- 
ing accreting matter into X-ray luminosity is ap- 
proximately 0.1 for neutron stars and black holes 



4 



( Lamers et all 19761) . The observed X-ray flux will 
also depend on the system distance, d, and the H I 
column density, nH. While the Lamers et al. model 
is useful for estimating the integrated X-ray flux 
produced by wind accretion, we emphasize that 
it is an approximation that does not include any 
spectral information about the emission. There- 
fore the predicted X-ray emission should be viewed 
with some caution. 

We could not measure the mass loss rate, M, 



from our optical spectra IjMcSwain et al.l I2007J ) 



since we did not observe the Ha line. Instead 
we obtained pu blished values of the observed 
M f or each star. | Chlebowski fc Garmanvi (|199lh 
and iHowarth fc Primal (|1989l ) found values for 
HD 14633 of log M= -6.81 M yr" 1 and 
—6.9 M Q yr -1 , respectively. We adopt the 
lower value for ou r mode l of the wind accretion. 
Howarth fc Prinia ( 19891 ) also measured log M= 
-6.5 M Q yr" 1 for HD 15137. 

We measured their terminal wind velocities, 
Voo, using flux-calibrated UV spectra obtained 
from the International Ultraviolet Exp l orer fI VE) 



archives. iPrinja. Barlow, fc Howarth! (|199(J ) de- 



scribe how to use either the time-variable nar- 
row absorption components of UV lines or the 
more steady-state, saturated line edge to mea- 
sure Vac, and we chose this latter technique due 
to the low S/N of the IV E spectra. For HD 
14633, we co-added the available high dispersion, 
Short Wavelength Prime Camera (SWP) spectra 
and looked for saturation in the N V A1238.821, 
Si IV A1393.76, and C IV A1548.202 lines. V x 
corresponds to the blue edge of the saturated line 
profile. HD 15137 has only one high dispersion 
SWP spectrum available, so we binned this by 5 
pixels to improve the S/N, reducing the precision 
of our measurement. The resulting values of Voo 
are listed in Ta ble [3 and they agree well with the 
values found by lHowarth et al. ( 1997 ). Knowledge 
of Voo enables us to characterize the stellar wind 
in the vicinity of the star according to the wind 
velocity law 



V(r) = Voo 1 - 



(5) 



with = 0.8 for weak s tellar winds (jHowarth fc Prinia 
19891 : IPuls et~allll996h . 

The remaining input parameters for the wind 
accretion model come from our orbital solutions, 



SEP fits, and stellar distances (jMcSwain et al 
2007). We used a grid of values for Mi with 
the minimum AI2 allowed by the mass functions 
(i = 90°). For these eccentric binary orbits, we 
expect the accretion rate to be greatest just af- 
ter periastron, where the wind density is higher 
and its velocity slower relative to the companion. 
Therefore we calculated the time-averaged X-ray 
flux, Fx, as well as the expected maximum and 
minimum flux during the orbit as a function of the 
primary mass. The results are illustrated in the 
bottom panels of Figures [1] and [5] We note that a 
higher M 2 increases the predicted X-ray flux since 
the gravitational force accreting the wind is larger. 

Normal O stars are expected to have an intrin- 
sic X-ray flux due to shocks in their stellar winds 
( Sana et al. 20061 and references therein). The 
winds are driven by the absorption and reemis- 
sion of UV photons in the spectral lines, and small- 
scale in stabilities in the wind will quickly grow into 
shocks (jLamers fc Cassinellill 19991 ). Magnetic con- 
finement of the stellar winds and colliding winds in 
close, massive bi naries may also co ntribute to the 
X-ray e mission dSana et al.l 120061 and references 
therein). Sana et al.l ( 2006F found that, with the 
exception of two colliding wind binaries, all of the 
binary and single O-type stars in the open clus- 
ter NGC 6231 display a tight correlation between 
their bolometric luminosities, L^i, and their X- 
ray luminosities, Ex, in the 0.5 — 10.0 keV band. 
Thus the intrinsic X-ray emission of an O star can 
be predicted from the relation, 

logix -logL boi = -6.912 ±0.153 (6) 



(ISana et al.ll2006j ). Usin g L^oi from the calibration 
of iMartins et all (|2005l ). we find that HD 14633 
and HD 15137 each have an intrinsic Ex of 2 x 10 31 
to 3 x 10 31 erg s " 1 . With the adopte d distances 



given in Table fTI ( McSwain et al.ll2007h . this corre- 
sponds to an unabsorbed flux between 3 x 10~ 14 
and 7 x 10~ 14 erg cm~ 2 s , a negligible contribu- 
tion to the predicted flux from stellar wind accre- 
tion. 

4.2. Observed Upper Limits 

Since these stars have not been detected as 
X-ray sources, we estimate an upper limit for 
their unabsorbed X-ray fluxes using other detected 
sources from the ROSAT All-Sky Survey and 
pointed ROSAT observations within a 30' radius 
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of each star (I White. Giommi. fc Angelinil 2000; 
Voges et alJl2000h . We expect that wind accretion 



will generate X-rays having a power law spectrum 
typical of MXRBs, and we used the photon in- 
dex r = 2.0 from the micro quasar LS 5039 (an- 



other wind-accreting MXRB. iBosch-Ramon et al 



2005). MXRBs in a lo w/hard emissio n state gen- 
erally have r ~ 1.6 ( Bellonil 2001 ). which wc 
also consider. The Galactic H I column density, 
nH, for each target was calculated using NASA's 
HEASARC toofj Finally, we used the HEASARC 
WebPIMMS calculatoJl to convert the upper limit 
of the X-ray count rate to the expected unab- 
sorbed flux in the 0.5—10 keV range. These details 
and the resulting upper limits for Fx are listed in 
Table|U We compare these upper limits to the pre- 
dicted Fx in the bottom panels of Figures Q] and 
[2j In each case, the observed upper limit is higher 
than the intrinsic X-ray flux from the O-star, con- 
sistent with the non-detection of each star. 

The wind accretion model predicts X-ray fluxes 
of HD 14633 and HD 15137 that are at least an or- 
der of magnitude greater than the observed upper 
limits for the most likely O star masses. However, 
it is worth noting that the wind-accreting micro- 
quasar LS 5039 exhibits order-of-magnitude vari- 
ations in both its X-ray flux and mass-loss rate s 



(|Bosch-Ramon et all Ex)! iMcSwain et~al"1l2004h . 



and othe r MXRBs may have flux variations of or- 
der 100 (|Bellonil[200lh . If the ROSAT observa- 
tions occured during times of low stellar mass loss, 
it is possible that any X-ray emission from a neu- 
tron star might be below the detection limit. O 
star winds may also atten uate any X-ray em ission 
in low energy bands, and ISana et alJ (|2004l ) find 
order-of-magnitude variations in the X-ray emis- 
sion from the colliding wind binary HD 152248 
due to variations in the absorption with orbital 
motion. If HD 14633 and HD 15137 are high incli- 
nation systems, they may experience similar vari- 
ations. 

4.3. Quiescent Neutron Stars 

Could these binaries contain neutron stars in 
a quiescent state? The conditions for wind ac- 



lr The HEASARC nH calculator is available online at 
|http: / /heasarc. nasa.gov/cgi-bin/ Tools / w3 nh/w3nh.pl| 

2 The WebPIMMS calculator is available online at 
http: / /heasarc. nasa.gov/Tools/w3pimms. html 



cretion onto a neutron star depend strongly upon 
the spin rate and the magnetic field, an d these are 
discussed in detail by iLipunovl (|l992h . A young 
neutron star will not accrete significant amounts 
of material because its fast rotation and large mag- 
netic field sweep material out of the system at a 
distance larger than the accretion capture radius, 
R a - Systems in this ejector regime are often ob- 
served as pulsars, and they are expected to remain 
in this state for 10 5 — 10 6 years. The spinning mag- 
netic field will dissipate rotational energy, causing 
the neutron star to spin down with time. The neu- 
tron star passes into the propellor regime, where 
the magnetospheric radius, R ma g, is smaller than 
R a but remains larger than the corotation radius, 
R c . The rapid rotation continues to inhibit accre- 
tion onto the neutron star. Eventually the com- 
pact star spins down enough to allow R c to expand 
outward beyond R m ag, trapping material in the 
vicinity of the neutron star and allowing accretion 
to begin. 

We can investigate whether the proposed neu- 
tron stars should be in the accretion re gime us- 
i ng th ese characteristic distances given bv lLipunov 
(|l992h . Assuming a nominal Ma = 1.4 M Q for the 



neutron star and using the relative wind/orbital 
velocity V re i ~ 1500 km s _1 for both systems, 
their accretion capture radius is 



Ra = 



2GM 2 
~V T 7 



2 * 10 10 cm. 



(T) 



Using the rotational frequency uj w 1 s _1 (dis- 
cussed in Section 3), we can estimate the corota- 
tion radius as 



Rr,= 



( GM 2 

\ LJ 2 



1/3 



6 * 10 8 cm. 



(8) 



The magnetospheric radius, R m ag, can be deter- 
mined by equating the ram pressure of the accret- 
ing gas to the magnetic field pressure. While we 
cannot measure the magnetic moment, n, the ob- 
served range among other neutron stars is 10 28 to 
10 31 ' 5 G cm 3 , with the mo st typical value being 
10 30 G cm 3 (jLipunovlll99"2h . Thus we find 



Rn 



2S a y/2GM 2 



2/7 



3 * 10 9 cm (9) 



for both binaries. Unrestrained accretion onto the 
neutron star requires that R mag < R a and R mag < 
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R c , which may be possible in these systems given 
the large uncertainties in R c and R mag - 

HD 14633 and HD 15137 could contain neutron 
stars in a low accretion state, and new X-ray ob- 
servations should conclusively show whether these 
systems harbor compact companions. Even qui- 
escent neutron stars are expected to produce an 
X-ray luminosity of 10 30 to 10 34 erg s - 1 with a 



thermal or nonthermal spectrum (jLipunov 1992), 
somewhat comparable in magnitude to the ex- 
pected intrinsic luminosity of an O star but pre- 
sumably distinguishable by the spectral energy 
distribution. Therefore any detection of X-rays 
with some spectral resolution should be able to 
distinguish between the soft spectrum of stellar 
wind emission, the presence of a quiescent neu- 
tron star, or emission produced by wind accretion 
with a power law spectrum. 

5. The Dynamical Interaction Scenario 

Both HD 14633 and HD 15137 contain a sec- 
ondary star with a minimum mass M% = 1 — 
2 Mq, so a neutron star companion seems plausi- 
ble. However, we do not detect pulsars in either 
system, and neither HD 14633 nor HD 15137 ex- 
hibits evidence of X-rays produced by accretion 
onto a compact companion. Although we have dis- 
cussed several reasons to account for the lack of ra- 
dio or X-ray detections, we must explore the pos- 
sibility that neutron stars are not present in these 
systems. We consider here whether these runaway 
SBls may have been ejected from NGC 654 by 
close dynamical encounters in the open cluster. 

In our earlier stu dy of these runaway stars 
( Boyaiian et al.l l2005h . we thought the dynam- 
ical interaction scenario unlikely based on the 
low mass ratios, q < 0.3, of HD 14633 and 
HD 15137. There are no known O stars in 
close, spectroscopic binaries that contain main 
sequence or pre-main seq uence companions late r 
than mid B-type stars ( Pourbaix et al.l l200-4h . 
and binary O-type stars usually form in open 
clusters with q ~ 1 ( Clarke fc Pringle! Il992l 
and references therein). Therefore it is unlikely 
that these binaries were born with low q and 
then ejected. However, the current mass ra- 
tios can be explained by dynamical interaction 
scenerios that exchange binary partne r s (e. £ 
Gualandris. Portegies Zwart. fc Eggletonl 2004 



The pairs may have formed in separate exchanges 
involving a wide O-type pair and a close B-type 
pair with a comparable binding energy. Such in- 
teractions can also explain the high e of the result- 
ing binaries. They would retain this eccentricity 
for more than 10 10 years since the low mass ra- 
tio significantly increases the timescale of tidal 
circularization of the orbit (Hil ditchll200l . 

To test the plausibility of this dynamical ori- 
gin, we simulated the ionization of a wide O-type 
binary in three- and four-body scattering interac- 
tions using the Starlab package developed by P. 
Hut, S. McMillan, J. Makino, and S. Portegies 
Zwart. In each simulation, we assumed a 20 Mq + 
20 Mq pair in a circular orbit with a semi-major 
axis of 1 AU. The radius of each O star was set 
to 8 Rq to represent a typical unevolved pair of 
late O-type stars. In our three-body simulations, 
we used a projectile consisting of a single 5 Mq 
star with radius 4 Rq and velocity of 140 km s -1 . 
The B-type projectile with this velocity at infinity 
would approach the O star binary with a kinetic 
energy, K, of only 0.13 times the gravitational po- 
tential energy of the binary, making an ionization 
extremely unlikely. In 50,000 three-body simula- 
tions in which we allowed the impact parameter of 
the projectile, r, to vary from to 5 AU, the target 
binary became unbound in fewer than 1% of our 
trials. For r > 1.25 AU, the chance of ionization 
decreased to nearly zero. Increasing K by factors 
of 2 and 3 (to achieve incoming velocities of 200 
km s _1 and 240 km s" 1 , respectively) reduced the 
number of ionizations further. 

A four-body interaction is far more likely to 
separate an O star binary. For these simulations, 
we used a target pair identical to that described 
above, and the projectile consisted of a 5 Mq + 
5 Mq pair of B-type stars in a circular orbit with 
a semi-major axis of 0.12 AU. We performed 8200 
simulations using the same range in K as above, 
thereby allowing incoming velocities of 100 — 170 
km s _1 for the projectile binary and allowing r to 
vary from to 5 AU. The O binary was ionized 
in 11% of our trials. We observe no strong depen- 
dence on the ionization rate for r < 5 AU with the 
lowest K, although higher K projectiles produce 
fewer ionizations when r > 3 AU (they tend to 
fly past the target binary, preserving the original 
configurations). 

In most of our four-body trials, Starlab pre- 
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dieted collisions between the targets and/or pro- 
jectiles, regardless of whether the target binary 
became separated. Without a robust hydrody- 
namical simulation, it is difficult to accurately pre- 
dict how many four-body interactions would re- 
sult in a close, bound pair in a binary-binary ex- 
change rather than a stellar collision. Neverthe- 
less, our simulations did predict a small number 
of exchanges (10) in the four-body scattering sce- 
nario. 

6. Conclusions 

Certainly it is plausible that HD 14633 and HD 
15137 were both ejected in multi-body encounters 
in the dense open cluster NGC 654. Simulations of 
a typical four-body interaction suggest that 11% 
of encounters between low mass (5 Mq + 5 Mq) 
binaries will ionize an O star pair (20 Mq + 20 
Mq). However, this rate is too low for such a dy- 
namical interaction to be probable, and we can- 
not conclude that these binaries were ejected from 
NGC 654 through such an interaction. 

While we do not presently detect a neutron star 
in either system, we cannot rule out that these 
systems were ejected by supernovae in the close 
binary systems. We recommend performing fur- 
ther X-ray observations deep enough to detect the 
intrinsic emission from the stellar wind shocks in 
these O stars. Any additional emission from a qui- 
escent neutron star or from active accretion onto 
its surface should be comparable to or stronger 
than the intrinsic wind emission, thereby produc- 
ing a measurable X-ray excess. Such observations 
are necessary to determine whether these systems 
contain neutron stars and thus determine conclu- 
sively which mechanism is responsible for produc- 
ing these runaway binary systems. 
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This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 
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Fig. 1.— (a) The mass diagram of HD 14633 
is plotted for a range of inclination angles (dot- 
ted lines). The most prob able values of M2 from 
iMazeh fc Goldberd (Il992h (solid l ine) a nd the ex- 



pected Mi from Martins et al.l ( 20051 ) (vertical 



dashed line) are also shown, (b) Predicted time- 
averaged X-ray flux (solid line) and the predicted 
range in flux over the eccentric orbit (gray) are 
shown for the lowest possible M 2 . The observed 
upper limits for Fx are shown for comparison (dot- 
ted lines). 
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Fig. 2. — The mass diagram and predicted X-ray 
flux are plotted for HD 15137 in the same format 
as Figure Q] 
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Table 1 

Summary of Results from McSwain et al. 2007 





P 




/(to) 


a\ sin i 


d 






Star 


(d) 


e 


(M e ) 


(Re) 


(pc) 


(K) 


logs 


HD 14633 


15.4082 


0.700 


0.0040 


4.13 


2040 


35100 


3.95 


HD 15137 


30.35 


0.48 


0.004 


6.7 


2420 


29700 


3.50 



Table 2 

Summary of Green Bank Telescope Observations 



Star 


MJD - 
Receiver 2450000 


Orbital 
Phase 


Exposure Su m 
time (s) (A*Jy) 


-Mim 

(mJy kpc 2 ) 


HD 14633 


S-band 3909.4194 
820 MHz 3909.5085 


0.381 
0.387 


6300 15 
6300 100 


0.062 
0.42 


HD 15137 


S-band 3897.4847 
820 MHz 3897.5753 


0.482 
0.485 


6300 15 
6300 65 


0.088 
0.38 


Table 3 

Terminal Velocity Measurements 


Star 


NVA1238.821 Si IV A1393.76 

Voo (km s- 1 ) Voo (km s" 1 ) 


C IV A1548.202 

Voo (km s- 1 ) 


Adopted Voo 
(km s _1 ) 


HD 14633 
HD 15137 


1677 
1660 


1710 


1690 


1677 ± 10 
1690 ±25 



Table 4 
Observed X-ray Flux Limits 







HD 14633 


HD 15137 


Value of nH (cm 2 ) 




6.76 x 10 20 


2.22 x 10 21 


Upper limit of count rate (cts s _1 ) 




0.00587 


0.0129 


Energy band of ROSAT /PSPC observation 


(keV) 


0.24-2.0 


0.1 - 2.4 


Assumed T 




2.0 


2.0 


Unabsorbed Fx, 0.5 — 10 keV band (erg cm 


" 2 s" 1 ) 


1.60 x 10~ 13 


5.62 x 10~ 13 


Assumed T 




1.6 


1.6 


Unabsorbed Fx, 0.5 — 10 keV band (erg cm 


" 2 S" 1 ) 


2.39 x 10" 13 


7.85 x 10~ 13 
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